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I. INTRODUCTION 

An important advance in the investigation of quantum 
fluids was recently achieved with the experimental obser- 
vation of high quantum degeneracy and ofF-diagonal long- 
range coherence, in a gas of exciton- polaritons in a two- 
dimensional semiconductor microcavitji^. High quan- 
tum degeneracy has been also observed in long-living 
polariton systems close to thermal equilibrium^. Bose- 
Einstein condensation (BEC) is the most natural way 
of describing these findings. However, due to the pecu- 
liarity of the polariton system, in particular the finite 
polariton lifetime, the intrinsic 2-D nature and the pres- 
ence of interface disorder, the existing theoretical frame- 
works rather interpret the phenomenon in strict analogy 
either with laser physics^i^ or with the BCS transition 
of Fermi particles^i^ii^. In particular, the problem of the 
polariton kinetics and of the non-equilibrium effects, due 
to the finite polariton lifetime and the relaxation bot- 
tleneck, have been investigated in many worka^'^ d'-''^^ . 
From these studies, the main effect of non-equilibrium 
seems to be a significant depletion of the condensate - 
with corresponding loss of long-range coherence — and 
the appearance of a diffusive excitation spectrum at low 
moment ai^iii. 

In spite of the high relevance of these theoretical de- 
scriptions, two basic questions remain still unanswered. 
Are the experimental findings correctly interpreted in 
terms of a quantum field theory of interacting bosons? 
Could the achievement of polariton BEC give new in- 
sights into the fundamental physics of interacting Bose 
systems? 

In this work we tackle these two questions, by show- 
ing that polaritons can be modeled borrowing from the 
theory of interacting Bose particles. In particular, we 
describe self-consistently the linear exciton-photon cou- 
pling and the exciton-nonlinearities, by generalizing the 
Hartree-Fock-Popov (HFP) description of BEC to the 
case of two coupled Bose fields at thermal equilibrium. 
In this way, we compute the density-dependent energy 
shifts and the phase diagram and we find a very good 



agreement with the recent experimental findings. Then, 
we apply the present theory to derive the full set of equa- 
tions describing the density-density response of the po- 
lariton condensate to an external perturbation. Focusing 
on the photon density response, which is directly related 
to photoluminescence, we predict different response of 
the collective modes to optical (light) or mechanical (co- 
herent phonons) perturbations. Since this behavior is 
driven by the presence of a coherent exciton field, we 
suggest that an experiment investigating these features 
could possibly solve the tricky problem of assessing the 
nature of the polariton condensate. 

II. THEORY 

The physics of the polariton system is basically that 
of two linearly coupled oscillators, the exciton and the 
cavity photon field n^'^d^ . Considering the limit of low 
density, the exciton field can be treated as a Bose field, 
subject to two kinds of interactions, the mutual exciton- 
exciton interaction and the effective exciton-photon in- 
teraction, originating from the saturation of the exciton 
oscillator strengthi^. Therefore, to describe polariton 
BEC we extend to the case of two coupled interacting 
Bose fields the formalism adopted in describing the BEC 
of a single Bose fiel d^^d^ . 

We express the exciton and photon field operators in 
the Heisenberg representation via the notation 

^,(r,i) = ^E^^''"^k(i), (1) 

and 

^,(r,t) = -^^e''^'-Ck(t), (2) 

where A is the system area, while 6k and Ck are inde- 
pendent Bose operators ([6k, cj^/] — 0). Notice that in 
this work we assume scalar exciton and photon fields. 
However, the theory can be generalized to include their 
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vector nature, accounting for light polarization and exci- 
ton spin^^. We adopt a finite system area A in order to 
model the effect of confinement, due to both the intrin- 
sic disorder— and to the finite size of the excitation 
spo t^°i^^ . While in 2-D, in the thermodynamic limit, 
the occurrence of BEC would be prevented by the diver- 
gence of thermal fluctuations^^, the finite size modifies 
the density of states, resulting in a finite amount of ther- 
mal fluctuations^^. The dependence of the results on A 
is discussed in Section Hill 

The exciton-photon Hamiltonian, including the exci- 
ton non-linearities, reads 

H = Hq + Hr + H^ + Hs. (3) 

The non-interacting exciton-photon Hamiltonian is 

^0 = E (^£^k^k + 44ck) , (4) 

k 

= h^k^ /2mx is the free exciton energy dispersion, mx 
is the exciton effective mass, = e^yjl + [k/kzY is the 
free photon dispersion, efj = h{c/nc)kz, c is the velocity 
of light, ric is the refractive index, kz = tt/Lc, and Lc is 
the cavity length. The term 

Hii = nn„"^{blcu + h.c.) (5) 

k 

describes the linear exciton-photon coupling. The term 
= ^ E ^'.(k,k',q)6^q&t,_qSkik (6) 

k,k',q 

is the effective exciton-exciton scattering Hamiltonian, 
modeling both Coulomb interaction and the non-linearity 
due to the Pauli exclusion principle for electrons and 
holes forming the exciton. The remaining term 

^■^ = i E ^^s(k,k',q)(ct+q6t,_q6kik + h.c.) (7) 

k,k',q 

models the effect of Pauli exclusion on the exciton os- 
cillator strengtbiS., that is reduced for increasing exciton 
densitySi. In this work we account for the full momentum 
dependence of t;2;(k, k', q) and Ws(k, k', q}^i^. In par- 
ticular, these potentials vanish at large momentum, thus 
preventing the ultraviolet divergence typical of a contact 
potentialii, without introducing an arbitrary cutoff. 



A. Bogoliubov ansatz 

To describe the condensed system, we extend the Bo- 
goliubov ansatsii to both the exciton and photon Bose 
fields, 

^x{c} (r, t) = ^x(c) (r, t) + i'x(c) (r, t) , (8) 



i.e. the total field is expressed as the sum of a clas- 
sical symmetry-breaking term ^x{c) for the condensate 
wave function, and of a quantum fluctuation field ipx{c)- 
The Bogoliubov ansatz imposes to introduce anomalous 
propagators for the excited particles, describing processes 
where a pair of particles is scattered inwards or outwards 
the condensate^i^. The resulting 16 thermal propaga- 
tors in the matrix form (in the energy-momentum repre- 
sentation, assuming a spatially uniform system) are 

G(k, ^u„)=( J' J' '""J ) , (9) 

where the elements of each 2x2 matrix block are (j, I = 
1,2; X,£.= x,c) 

gffiKicu^) = - [^dre^--^d{ (k, r) (k,0)^).,;3 , 
Jo 

(10) 

hcun = 27rn//3, n = 0,±1,... are the Matsubara energies 
for bosons, /3 = l/ksT and the symbol {...)t,i3 indicates 
the thermal average of the time ordered product. Here, 
to represent the exciton and the photon fields, we adopt 
the compact notation 6^(k) = (^^(k), d|(k) = 0|(-k) 

and 0:r(k) = 6k, Oc(k) = Ck. Correspondingly, the gen- 
eralized one-particle density 

(11) 

with X, ^ = x, c, is separated into the contribution of the 
condensate rv^^ = $J<&^ and of the excited particles 

"^^ = E"^^-E<Ox(k)0|(k)). (12) 

This latter quantity represents the excited-state density 
matrix, expressed in the exciton-photon basis, and it is 
directly related to the corresponding normal propagator 
via the well known relationi^ 

= - / ^Im{(gf«)'-*(k,c.)}nB(c^) , (13) 

where the retarded Green's function 

(g?«)'-^*(k, oj) = ,g?«(k, iujn^oj + zO+) (14) 

is the analytical continuation to the real axis of the 
imaginary-frequency Green's functioni^. 

B. Condensate wave function 

Within the Popov approximation, for a uniform sys- 
tem, the two coupled equations for the condensate am- 
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plitudcs are 

ih^., = [eg-2Rc{w,(0,0)ng= + 2^'w,(k,0)nn 

k 

+ («,(0, 0)nS- + 2 ^' 5,(k, 0)nr)]<i>x 

k 

k 

-[vs{0, 0)ng^ + 2 ^' «,,(k, 0)niJ^]$, (15) 



k,0 



where = Ek^o ^nd 



(k, q) = 2 Kc'') °) + ^^{*) ~ ^)] ■ 

(16) 

We assume that both the condensate fields evolve with 
the same characteristic frequency E/h, i.e. 



(17) 



By replacing this evolution into Eq. (jlSp . we obtain a 
generalized set of two coupled time-independent Gross- 
Pitaevskii equations, which can be formally written in 
the matrix form 



p ( -^0 \ _ fGP f ^0 



(18) 



where 



Lff ^ e^- 2Rc{«,(0, OX + 2 ^ ' f;,(k, 0)7i£^(}19) 

k 

+ («,(0,0)ng^ + 2^'«,(k,0)nr) 

k 

if/ = aOfl^-^f;,(k,0)n£- 

k 



-GP 
^22 — "=0 ' 



(20) 



we have defined the normalized Hopfield coefficients of 
the condensate state as = -'^o'i' and <I>c = Co$, sat- 
isfying |Xop -I- |CoP = 1, and rig = is the actual 
density of the polariton condensate. The two solutions 
E = E^P^^'P'^ of Eq. (UHl) define the lower and upper po- 
lariton condensate modes 



(21) 



The condensate energy is given by the lower energy solu- 
tion E^Q , which corresponds to the minimal energy of the 
polariton states. In the present U(l) symmetry-breaking 
approach, we can identify the condensate energy with the 
chemical potential of the polariton system, i.e. E^^ = /i. 
The grand-canonical thermal average has to be taken ac- 
cordingly. 



C. Beliaev equations 

In analogy with the standard field theory for a single 
Bose fielcU^, the 4x4 matrix propagator (^(k, iuin) obeys 
the Dyson-Belaev equation 

G (k, iiOr,) - G" (k, iuJn) [1 + S (k, iUn) G (k, iW„)] , 

(22) 

where we have introduced the matrix of the non- 
interacting propagators 

G° ^ {g° (k,zc.„)}^« ^ 6Mi-y^^n^^k^+iA-' (23) 
and the 4x4 self-energy matrix 

■ S--(k,ic^„) E"^(k,zw„) ■ 



I](k, ILOri 



S'=-(k,zc^„) E^^(k,zw„) 



(24) 



Within the HFP limit, the self-energy elements are inde- 
pendent of frequency and read 

S^;(k) = 2^ [^^.(k,q)n^- - t;.(k,q) « -!-<)] , 
q 

S^2"(k) = (S^f )* = «.(k, ml - 2^^s(k, 0)$,$,, 
S-(k) = S]2-(k) = hilR - 2^«,(k,q)n^-^ , 



= (sff(k))* = -vs^l. 



(25) 



while S^"f (k) = SJf(-k) and Sj=^(k) = 0. 

The solutions of Eq. ([22]) can be written analytically 
in terms of the self-energy elements and the unperturbed 
propagators. For example we obtain 



ff?f(p) 



.gg(p) [l-gg(-p)jV£,(p)] 
|l-5g(p)7Vz,(p)|^-|<7g(p)iVB(p)r 



where p = k, iun, 

Nd{v) - STf(k)+gg(p)|E?^(k)p+gS(-p)|ET2^(k)p, 
A^b(p) = STf(k) + [.9S(p)+.9j)(-p)]Srr(k)Sf2^(k), 
and 

^ _ gg(-p)iV^(p) 

For each value of k, the analytic continuation of each 
Green's function g^^{]<i,z) shares the same four simple 

poles at z = ±e'^^^^\ i.e. 



E - 



x^(k)p , 



j=lp,up 



Ei{k) z + Ei(ky 



g-,^{k,z) = J2 



j=lp,up 



z-Ei{k) z + E^{k)* 



|C^(k)|2 , |C^-(k)|^ 



j=lp,up 



-E3{k) z + Ei{k)* 



_ ^ xi{kYCi{k) , xjikyciXk) 

5n(k,^j - 2. z-&{k) + z + Ei{kY ' 

j=lp,up 
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and so on^i The poles of the propagators represent the 
positive and negative Bogohubov-Behaev eigen-energies 
of the lower- and upper-polariton modes. The residual 
in each pole depends on the corresponding generalized 
Hopfield coefficients. 

We point out that the polariton excitation modes for 
a given k can be also obtained by directly diagonalizing 
the problem 




L 



HFB 




(28) 



with 



^HFP 



and e^^'^^ 



+ ^11 


\^xx 
^12 


^11 


^12 


\^xx 
^21 


^l^k + ^22 ) 


•^xc 
^21 


~^22 


^11 


^12 


^k-M 





\^xc 
~^21 










(29) 



^k'"'^'' ~ "^^^ components of the 4 eigen- 



vectors hj(fc) EE {Xu,Xy,Cu,Cy)j{k) (j = 1,...,4) are 
again the generalized Hopfield coefficients corresponding 
to the normal {Xu,Cu) and anomalous (Xi,,C„) compo- 
nents of the polariton field, in analogy with the one-field 
HFP theory. They obey the normalization relation 



1 



(30) 



a condition which guarantees that the operator destroy- 
ing the lower (upper) polariton excitation with wave vec- 
tor k, 



tI = X^(k)Sk + X^(k)6lk 

+ C^„(k)6k + C^(k)cLk: 



(31) 



j = Ip, up, obey Bose commutation rules. The lower (up- 
per) polariton one-particle operators are then defined 
by 



u^'(k)pk + v^i-kypl 



k ' 



(32) 



where the normal and anomalous polariton coefficients 
are given by 

v^ik)^[Xiik) + Ci{k)f\ 
W{k)\'~\v^{k)\^ 



1 . 



(33) 



The normal modes of excitation are thermally populated 
via the Bosc distribution 



1 



1 



(34) 



while the lower- and upper-polariton one-particle densi- 
ties are given by 



K 



1 

A 



(|^.^(k)p + K(k)p)iV^+K(k)|^ 



(35) 



The first and the second term of the sum represent the 
thermal and quantum fluctuations, respectively. There- 
fore, for a fixed total polariton one-particle density Up, 
the density of the polariton condensate is given by 



no = |$r = np 



kT^O 



(36) 



From no, the exciton and the photon condensed densities 
are finally obtained via Eq. ([20)l . 

Hence, for a given polariton density Up and temper- 
ature T, a self-consistent solution can be obtained by 
solving iteratively Eqs. (HH]), ^ and until 

convergence of the chemical potential fi and the density 
matrix n^j(k) is reached. From this self-consistent so- 
lution, we obtain the exciton and photon components of 
the condensate fraction as well as the spectrum of col- 
lective excitations and the one-particle populations. We 
point out that the self-consistent solution must be in- 
dependent on the initial condition used in Eq. (|22p and 
([TO]). Fast convergence of the iterative procedure in the 
numerical calculations is obtained by starting from the 
ideal gas solution, i.e. the solution obtained by neglecting 
the two-body interactions, and considering the resulting 
polariton states occupied accordingly to the Bose distri- 
bution. 



III. RESULTS AND DISCUSSION 

For calculations we adopt parameter modeling typical 
GaAs based microcavity samples^i^i^, in particular we 
assume a linear coupling strength hflu = 7 meV, corre- 
sponding to 12 embedded GaAs quantum wells, and the 
photon-exciton detuning S — — €q =3 mcV. Where 
not differently specified, we consider a system area A = 
1000 /im^ and a polariton temperature T = 10 K. For 
the interaction potentials Vx and Vs, we use momentum- 
dependent values following Rochat et a\^. 



A. Spectral and thermodynamic properties 



In Fig. [l]we show the energy- momentum dispersion of 
the collective excitations, ii?^ and ±£'^^, as obtained 
at the critical polariton density = 5 ^mT^ and far 
above the critical density, i.e. Up = 50 ^m~^. The 
curves correspond to the positive- and negative-weight 
resonances for the lower- and the upper-polariton. We 
notice that, for the largest value of rip, the polariton 
splitting decreases, due to both the exciton saturation, 
decreasing the effective exciton-photon coupling SJf , and 
the change of the exciton-photon detuning produced by 
the exciton blueshift, given by Y,ff. However this varia- 
tion is quantitatively small, suggesting that, at equilib- 
rium, the polariton structure should be robust even far 
above the condensation threshold. We also mention that, 
close to zero momentum, the dispersion of the lower po- 
lariton branch, above threshold, becomes linear, giving 
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FIG. 1: The dispersion of the normal modes of the system 
for polariton density Tip = 5/im~^ (solid) and rip = 50/im~^ 
(dashed) . The uncoupled photon (dotted) and exciton (dash- 
dotted) modes are also shown. 
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FIG. 2; Lower (solid) and upper (dashed) polariton energies 
at = vs polariton density rip. Inset: Double logarith- 
mic plot of the lower polariton energy. The thin dotted lines 
highlight the two different slopes below and above the density 
threshold. 



rise to phonon-like Bogolubov modes, as in the standard 
equilibrium single-field theorj*ii. The modification of the 
energy splitting between the lower and the upper polari- 
ton branch is accurately characterized in Fig. [21 where 
the energy shifts of the two polariton modes at k = are 
plotted as a function of the density. Exciton saturation 
and interactions result in a global blue-shift of the lower 
polariton and a red-shift of the upper polariton. The 
shifts are linear as a function of the density, but their 
slope varies close to threshold. As highlighted in the in- 
set, the slope of the lower polariton shift changes by a 
factor of two across the threshold, because the contribu- 
tion of the condensed populations (n°^, n°^) is one half 
the contribution of the thermal populations {fixx^ n^c), 
as seen in Eq. 

We now turn to the thermodynamic properties of the 
system. In Fig. ^a), we report the density-temperature 
BEC phase diagram, as computed for A ~ 1000 /xm^. 
The phase boundary in the calculations has been set by 



the occurrence of a finite fraction of polariton condensate 
larger than 1%. In the plot, a few values of the quan- 
tity |Xop along the phase boundary are indicated. This 
quantity represents the exciton amount in the polariton 
condensate. It decreases for increasing density, depend- 
ing on the exciton saturation and the change in detuning. 
For very large densities this quantity eventually vanishes, 
corresponding to the crossover to a photon-laser regime. 
However, for the studied GaAs model microcavity, the 
variation of the exciton amount in the condensate field 
remains very small up to densities far above the exper- 
imentally estimated polariton density^i^. This is ba- 
sically due to the positive cavity detuning, sufficiently 
large to be robust to the exciton nonlinear energy shift. 
On the other hand, due to the same feature and to the 
flat energy dispersion of the cxciton-likc lower polariton 
states, a large population can be accommodated in the 
excited state when the temperature exceeds 25 K, thus 
dramatically increasing the BEC transition density, even- 
tually leading to the direct occurrence of photon-lasing. 
In particular, for this system, we predict that equilib- 
rium polariton BEC is impossible for temperatures larger 
than 30 K. In Fig. [3]Jb), we show a detail of the low-T 
region of the phase diagram, computed for two differ- 
ent system areas A = 100 /im^ and A = 1000 /im^. In 
a homogeneous two-dimensional system, in the limit of 
infinite size, a true condensate cannot exist due to the 
divergence of low-energy thermal fluctuations. The tran- 
sition to a superfluid state is instead expected, giving 
rise to the Berezinski-Kostcrlitz-Thouless crossover with 
spontaneous unbinding of vortices. The divergence of the 
condensate fluctuations has however a logarithmic depen- 
dence on the system size. Fig. [SJb) shows this behaviour 
as a slow increase of the critical density for increasing 
A. Quantitatively, the critical density varies by no more 
than a factor 2 at T = 1 K, for the two considered values 
of the system area. This difference becomes even smaller 
for larger temperatures. The predicted dependence on 
the system size could be experimentally verified only in 
samples with improved interface quality and manifesting 
thermalization at very low polariton temperature^^. 



IV. LINEAR RESPONSE IN THE HEP LIMIT 

Within the present theory, it is easy to compute the 
density fluctuation of the exciton or the photon field pro- 
duced by a perturbation acting on either of them. This 
quantity is particularly interesting because as we sug- 
gest later, it might be used to study the nature of the 
polariton condensate via fully optical or mechanical per- 
turbations. 

We consider the time dependent perturbation 



(37) 



driven by the external potential Vext{^,t) affecting the 
field 5'^(r, t) (x = x,c). This perturbation results in a 
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FIG. 3: (a) Phase diagram of the polariton condensation, 
computed for A — 1000 fitn^. The exciton fraction in the 
condensate |^op, along the phase boundary, is indicated in 
boxes, (b) Detail of the low-T region where we also display 
the transition boundary computed for A = 100 /xm^. 



density fluctuation 

fc««(r, t) EE n«(r, t) - n«(r, t) , (38) 

= x, c) around the equilibrium value ri,^|(r, i). 
Within the hnear response limit, this fluctuation is given 

6n^^ir,t) ='-j dt' J dr'Qx(r',t';r,t)K.t(r',0, 

(39) 

where 

C^^ir', t'; r,t) = { [h^^{r', i'), n«(r, t)] ) . (40) 

For a spatially uniform system in steady state, the den- 
sity commutator Eq. ([^0]) only depends on r — r' and 
t — t' . Then, by Fourier transforming Eq. (|39p . we get 
the expression 



where 



i^yk,.).-!/^^^^;^ (42) 

n J 2-!^ (jj A- «0+ 

is the retarded density-density correlational, and 
Cjj^(k, w) is the Fourier transform of C^^(r, i). 

Using the framework developed in Section [TTl we write 
the real-time density operators fi^^ (r, t) in the polariton 



basis, via Eqs. (IT][2]) and by inverting the transformation 

Eq. (inD 

j=lp,up 

(43) 

Here, the factors n^(k) and T^(k) define the compo- 
nents of the exciton (photon) field on the forward and 
backward propagating lower and upper polariton cigen- 
modes. Within the HFP limit, the density commutator 
consists of three contributiona^S^ 

Qx(r,t) = °C|^(r,i) +° Cgir,t) + C^^{r,t) (44) 
= [^m^{t){[4>dO, 0), (r, t)] ) - h.c. 



($^(0)$*(t)([v;«(O,0),^x(r,t) 



h.c 



+ { [^{0, 0)^5(0, 0), ^j^(r, t)^^(r, i)] ) , 

the first two terms arising from the presence of the con- 
densate fields. Correspondingly, the retarded density- 
density correlation can be written as the sum of three 
terms 

Dl^{k,u) = ^Dl^{k,Lo)+''Dg{k,u) + b^^{k,oj) . (45) 

The first two terms survive only in the presence of a con- 
densate while the third one only depends on the thermal 
population (however it is affected by the modification of 
the one-particle spectrum induced by the condensate). In 
detail, the first term describes the excitation of particles 
out of the condensate and it is given by 



with 



, (46) 



^|,(k,t 



$g$*nf (k)nj^(k) + $|$^Tf (k)TJ^(k) 
huj- eI+ iO+ 



(47) 

The second term describes the de-excitation of the ther- 
mal population into the condensate and it is given by 



j=lp,up 



, (48) 



(41) with 



^L(k,c.) 



$5$xnf (k)TJ,(k) + $*$*Tf (k)nJ,(k) 



hcj- El + iO+ 



(49) 

The third term describes the oscillations of the thermal 
population and it is given by 



D^AKu;)= [^^(k,^) 

j=lp,up 



T|,(k,-^) 



(50) 



7 




0.04 

E 

I 0.02 





-d" 
E 

' -0.02 



-20 



(a) 



1 , 1.5 

k„(^m-') 



-10 10 

fico (meV) 




1 , 1.5 



J 










(c) 











](d) 



20 -20 -10 10 
hcD (meV) 



0.04 . 



0.02 g^; 

X 
3 

(D 

-0.02 

B 



■0.04 



20 



FIG. 4: Imaginary part of the retarded density-density correlations D'^^{'k,Lj) (panel a and c) and D]^^(k,Lo) (b and d), as 
computed for T = 10 K and = 5 fim~^. In panels (a) and (b) we display the k — lo dependence of the correlation in grey 
tones. In panels (c) and (d) the same quantities are plotted as a function of the energy hco and at the wave vector k — 0.5 iim~^ . 
In panels (c)-(d), the dashed line represents the contribution Z)^^(k, tj) arising from the oscillations of the non condensate. 
Delta peaks arise from the oscillation of the condensate. At high energy (i.e. at the upper polariton energy), the photon-photon 
response corresponds to absorption at positive energy and gain at negative energy, while the photon-exciton response has the 
opposite behavior. 



with 



T|,(k,.) 



E 



(k,q) (iV^-iV^_k 



R 



and 



(51) 



F|^(k,q) = n|(q-k)ni*(q-k)n^*(q)n;^(q) 

+ n|(q- k)Ti*(q - k)4*(q)T!^(q) ,(52) 

i?|;^(k,q) = ^(q-k)ni*(q-k)T|(q)T^*(q) 

+ ^(q- k)Ti*(q- k)T|(q)n!;(q) .(53) 



We use these equations to study the fluctuation Sricc of 
the photon density, directly related to the photolumi- 
nescence measured in experiments, produced by an op- 
tical (affecting the photon field) or mechanical (affect- 
ing the exciton field) perturbation. Wc take, as ex- 
ternal potential, a plane wave with wave vector kext, 
dclta-pulscd in time, i.e. Vexti'Tjt) = Voe^^"'*''^S{t — to). 
In this case, from Eq. (|^T|) . we see that Sn^^{\i,ijj) = 
Vb-D^^(k, w)(5(k — kext), i-e. the response is diagonal in 
k and is simply proportional to the correlation . The 
imaginary part of describes the energy transfer to the 
system and thus it is the most relevant function. 



The resulting quantities — Im{D^^(k, w)} with % = c, x 
are shown in Fig. [H We display the results for the 
photon-photon (panels (a)-(c)) and photon-exciton (pan- 
els (b)-(d)) correlation in the condensed regime. Poles 



at negative energy are present for both quantities, due 
to the spectrum modification induced by condensa- 
tion. While the condensate contribution defines collec- 
tive modes with infinite lifetime (delta-peaks), the con- 
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tribution from the thermal population is responsible for 
a finite linewidth, as can be argued from the equations 
above. The surprising result of this analysis is however 
the unusual behavior manifested by the photon-exciton 
response at high energy. In this case, the collective mode 
corresponds to gain at positive energy and to absorption 
at negative energy. Conversely, for the photon-photon 
response, we expect to observe the opposite behavior. 
This feature is due to the phases of the Hopfield coeffi- 
cients. As shown by Eqs. (|47l49p . the opposite nature of 
the collective modes generated by optical or mechanical 
perturbations would be observed in experiments only if 
an exciton coherent field is present. Furthermore, the 
relative amplitude of the two responses is proportional to 
the fraction of coherence of each field, giving direct access 
to the amount of exciton condensate. In this respect, we 
point out that, although any kind of perturbation would 
affect simultaneously both the exciton and the photon 
field^ii, the geometry of the system can be chosen in such 
a way that the effect of the perturbation on one of the two 
fields be dominant (see for example the static perturba- 
tion affecting the exciton field used by Balili et al^) . We 
thus suggest that an ideal tool to observe these features, 
would be a pump and probe experiment where the probe 
be in turn optical or mechanical (for example produced 
by a coherent acoustic waves^S). 



ing the BEC of microcavity polaritons in very close anal- 
ogy with the BEC of a weakly interacting gasi^i^l. In par- 
ticular we treat simultaneously both the linear exciton- 
photon coupling and interactions. We account for the 
presence of a non condensed population as well. Within 
this description we are able to predict the modification 
of the spectrum and the thermodynamic properties for 
increasing density. Since the theory allows to describe 
simultaneously the properties of the polariton, the pho- 
ton and the exciton fields, it allows the understanding of 
typical optical measurements. In particular, our analy- 
sis supports the interpretation of the recent experimen- 
tal findingsii^i^ in terms of BEC of a trapped gas. We 
have applied the theory to compute the density-density 
response of a polariton gas to an external perturbation. 
This quantity can be characterized in pump and probe 
experiments realized with an optical or mechanical probe. 
In particular, we predict that the upper polariton energy 
collective modes generated in the two cases have a re- 
sponse of opposite sign. We suggest that the observation 
of this feature would be a proof of the presence of an 
exciton condensate and would answer the long standing 
question about the connection between polariton BEC 
and a laser phenomenon. 



CONCLUSIONS 



We generalized the HFP theory to the case of two cou- 
pled Bose fields at equilibrium. The theory allows modcl- 
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Shelykh et aZ.— Tmve recently studied the effects of po- 
larization and spin at T = 0, within the Gross-Pitaevskii 
limit restricted to the lower polariton field 
Here we write the general expression with complex energies 
E''P^^p'> (k). Within such a notation, the formulas can be in 
principle extended to include a phenomenological imagi- 
nary part to the energies, in order to account for the finite 
radiative lifetime of polaritons. 

We are only considering the response of the exciton and/or 
the photon density to a perturbation affecting one of the 
two species. For this reason we are not considering expres- 
sions involving the non diagonal terms of the density op- 
erator h^". 

In the HFP approximation, the coupling between the den- 
sity fluctuations of the condensate and the non condensate 
is neglected, because the HFP ground state is defined as 
the vacuum of quasiparticles. For a trapped gas of atoms, 
this approximation is found to be in qualitative agreement 
with experiments, although, close to Tc, deviations are 
reported. A better quantitative agreement has been ob- 
tained beyond the HFP limit, by means of Random Phase 
Approximation technique a^^'^^ . For our present purposes, 
however, the HFP limit is an adequate approximation. 



